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ABSTRACT 

The present work deals with the determination of effective elastic properties of fiber graded coniposite lamina, 
with varying fiber-volume fraction (FVF) along the thickness direction. In the available analysis, the effective elastic 
properties are evaluated based on continuous FVF and decoupled representative volume element (RVE) at every point 
along the thickness direction. A sub-volume is designed with varying FVF, according to sine function of thickness 
direction. In order to determine the effective graded material properties without any of these assumptions, a continuum 
micromechanics finite element model is derived. The RV is composed of several siib-volumes of different FVF and 
homogenization technique is applied over these sub-volumes, without decoupling these from the whole RV. The 
characteristics of effective elastic properti.es according to sine function are discussed. 
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INTRODUCTION 

In laminated composite structure design, delamination is the major drawback that limits the applications 
under high mechanical stresses and in thermal environment. In order to mitigate the delamination effect, the 
concept of functionally graded materials (FGMs) emerged, which is being widely used in designing advanced 
structures in various engineering applications. FGMs are heterogeneous materials in microscopic scale and 
attributed by continuous variation of material properties particularly in the thickness direction. The concept of 
FGM utilizing metal and ceramic as constituent materials was first proposed by Koizumi [1]. The FGM [1] is a 
particulate composite, in which, gradient of single composition is created along the thickness direction through 
exponential or power-law distribution of metal and ceramic particles. The properties of these FGMs are computed 
by any of these homogenous methods namely Mori-Tanaka, exponential distribution and power-law distribution. 
Martin and Leissa [2] first proposed the concept of graded fibre composite by varying the fibre volume fraction 
(FVF) along the thickness direction. 

Kubiak [3] considered variable FVF across the width of a thin-walled rectangular composite plate and 
investigated its dynamic response under the in-plane pulse loading. Batra and Jin [4] considered a graded 
composite plate of variable fiber-orientation angle along the thickness direction and investigated the effect of 
variation of fiber-orientation angle on the resonant frequencies of the plate. Cho and Rowlands [5] presented an 
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optimized local fiber-orientation angle, in order to reduce the stress concentration near the holes or notches of composite 
structure. Panda and Ray [6] considered laminated fiber-reinforced composite plates of variable fiber-orientation angle, 
along the thickness direction and analyzed their controlled nonlinear transient responses. Yas and Aragh [7] presented free 
vibration characteristics of cylindrical panels with continuous variations of FVF and fiber-orientation angle along the radial 
direction. Tahouneh et al.[8] analyzed free vibration characteristics of a thick annular continuously graded fiber-reinforced 
composite plate of variable FVF along the thickness direction. Aravinda et al [9] has evaluated the graded properties of 
composite lamina using a power law. 

A considerable number of work on graded fiber-reinforced composite structures has been reported in the 
literature. In most of the reported literature, the material properties are varied along the thickness direction by varying the 
FVF or by fiber orientation angle. In most of the available analyses on fiber-reinforced composites, the graded material 
properties are evaluated using standard micromechanical theories with an assumption of discrete representative volume 
element (RVE) at every point along the thickness direction. In the available literature, the properties corresponding to an 
RVE are computed by decoupling the RVE from the whole composite. The considerations like existence of RVE at every 
point and continuous variation of FVF cannot be physically achievable. Also, the discrete RVE will not ensure that, no 
effect on the graded material properties of the composite in macro scale. These assumptions in deriving graded material 
properties may lead to erroneous elastic responses of graded composite laminated structures. 

In the present work, microstructure of a graded fiber reinforced composite lamina with varying FVF along the 
thickness direction is considered. In order to obtain the graded material properties of the overall composite, 
micromechanics finite element analysis is carried out to predict the graded elastic properties of the composite lamina, 
without considering any of the assumptions reported in the earlier literatures. Presently, sub-volumes are stacked along the 
thickness direction. Each of the sub-volume has a continuous fibre with certain FVF. But, FVF discontinuously varies 
among the sub-volumes which follows the power law. Using micromechanics theory in conjunction with finite element, 
classical homogenization is applied on the sub-volumes without decoupling from the whole RV. Since the homogenization 
is applied on the coupled sub-volumes of different FVFs, the effect of their interactions is accounted. 


Graded Composite Representative Volume Element (RVE) 



Figure 1: Schematic Diagram of Representative 
Volume Element (RVE) 

The Figure.l shows the representative volume element of the graded lamina with different sub-volumes. 
Each sub-volume consists of a continuous fiber with certain FVF. The FVF at the middle sub-volume is maximum 
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compared to the FVF other sub-volumes which are stacked along +z direction and -z direction. The fibers are aligned along 
the x-axis and sub-volumes are stacked along the z-direction. The sub-volumes considered are square in cross-section with 

ac and the area of each fiber is a k f = a c yjv} . 


Finite Element Model of RVE 

In this section, a three-dimensional finite element model of RVE is presented. The strain vector at any point 
within the RVE can be expressed as, 

{£} = {£ £ £ Y Y Y } T , 

I J L xx yy zz / yz / xz / xy J ’ 

du dv r) w r)u r)v du r)w r)v r)w 

£ ,< = y . £ yv=T . £ zz = T ’ Yxy = 3— + T . ^« = T + T ’ ^"T + T (1) 

dx dy dz dy dx dz dx dz dy 


where, £ xx , £ and £ are the normal strains along x, y and z directions, respectively; J x _ and y are the 

transverse shear strains in the XZ and the yz -planes, respectively; y is the in-plane shear strain in the xy plane; 

ll(x, y, z), v(x, y, Z.), w(x, y, Z.) are the displacements at any point in the RVE along x, y and z directions, 
respectively. Similar to the strain vector (Eq. (1)), the stress vector at any point within the RVE can be expressed as, 

{(j} = {(T O <J T T T } T (2) 

1 J ^ xx yy zz yz xz xy J v ' 

where, <7 VV , <7 and <7 are the normal stresses along x, y and z directions, respectively; T x _ and T are the 

transverse shear stresses in the XZ and VZ -planes, respectively; T xy is the in-plane shear stress in the XV -plane. The 

displacements (u, v and w) at any point within the RVE can be expressed in the form of a displacement vector ({<:/}) as 
follows. 


form. 


{d} = {u V w} J (3) 

Introducing Eq. (3) in Eq. (1), the strain vector at any point within the RVE can be expressed in the following 

{£} = [L]{d} 


[L] = 


— 0 0 0 — — 


r)x 


dy 


dz dy 

± 0 — 
dz dx 


0 0 dz dy dx 0 


(4) 


constitutive relations for fibre and matrix phase materials within the RVE can be written as, 
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{^} = [C f/ ]{£-}, g= 1,2 (5) 

where, the superscript q denotes the quantities within the fiber or matrix phase volume according to its value as 1 
or 2, respectively; the matrix [C cj \ is the stiffness matrix for the fiber or matrix phase martial. The first variation of 
internal energy of the RVE can be expressed as, 

SU = j^\ [SeY{a q }dV t ( 6 ) 

q =i " 

where, 8 is an operator for first variation; V is the volume of fiber (ij = 1) or matrix (q = 2) phases. 

Introducing Eqs. (4) and (5) in Eq. (6), the first variation of the internal energy (SlJ ) can be written as, 

su = X f [Sd} T [L] T [C q ][L]{d}dV q (1) 

<7=1 V “ 

The volume of RVE is discretized by 27 nodded iso-parametric hexahedral elements. Accordingly, the 
displacement vector ({<:/}) for the i' h node of an element can be expressed as, 


K} = k M, W, } T , / = 1,2,3...,27 

(o) 

The displacement vector at any point within a typical element can be written as, 

M=WK) (9) 

where, [ /V] is the shape function matrix and [d } is the elemental nodal displacement vector. Using Eq. (8) in 
Eq. (9), the simplified expression for the first variation of internal energy ( 5U e ) of a typical element can be expressed as, 

SU-={Sd‘}\K-]{d‘}. [x'] = { v; ([A'] I [i] T [c 9 ][z.][M]), 7 v; (10) 


where, V ‘ is the elemental volume within fiber (q = 1) or matrix ( q = 2) phase volume. Upon assembling the 

elemental equations (Eq. 10) into the global space, the global expression of the first variation of internal energy can be 
obtained as, 

SU={SX} J [K]{X} (11) 

where [ K \ is the global stiffness matrix; {X} is the global nodal displacement vector. The internal energy of 
RVE (Eq. (11)) arises basically due to the application of homogeneous displacement boundary conditions over the 
boundary surfaces of the RVE. For a specified nodal displacement (say, X ; ), the first variation of the corresponding 

element of the global nodal displacement vector ({X}) is equal to zero (SX j = 0). Thus, the j th row of [X] is to be 
deleted and a column of [X] with the same index ( j ) (say, {/;}) is also to be removed for constituting the displacement 
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SU={Sxf([K r ]{X r }+(/,}*,) 


( 12 ) 


In Eq. (12), | K r ] and {X r } are the resulting stiffness matrix and nodal displacement vector after the imposition 
of specified nodal displacement. Equation (12) can also be written in generalized from when a number of nodal 
displacements ( N d ) are specified as follows. 


<St/={«■/([ AT r ]{X, }-{/>}) (13) 

Applying the principle of minimum potential energy i.e. SlJ = 0, the following governing equations of 
equilibrium can be obtained, 

[K r ]{X r } = {P} (14) 


For the applied displacement boundary conditions over the boundary surfaces of the RVE, the corresponding 
displacement field within the RVE can be obtained by the solution of Eq. (14). Utilizing this displacement field, the 
volume average strain vectors of RVE, sub-volumes, fiber phase and matrix phase can be computed by the following 
expressions, 



' (15) 

where, N* and N r/ - are the numbers of elements within the volumes of k th sub-volume and its fiber phase, 
respectively; N t is the total number of elements within the RVE; [ F t } is the strain vector at any point within the l th 
element; V} is the volume of \ th element. 


Table 1: Boundary Conditions for Determining Elements of {<? } 


Boundary Conditions 

Elements of {c } 

u 

w 

= 0, u 

—x 

= 0, w 

-z 

, =(4xC),vI =0, v| =0, 

+x v x c ' ’ 1 —y ’ \y 

+Z =0 

£ ~£°,£ = £ = f = r = f =0 

u 

w 

= 0, u 

—x 

= 0, w 

-z 

+x =0,v|_ y =0, v | y = ( E ° y xa c ), 

= 0 

z 

£ ~ £°,£= £ = f = f = f =0 

u 

w 

= 0, 

— X 

= 0, w 

-z 

u\ =0, v =0, v =0, 

\+x v—y \—y 

+ z =( £ °z xh c ) 

£ ~ £°,£ =£ =f = f = f =0 

z z ’ x y i yz i xz i xy 
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Table 2: Effective Elastic Coefficients of an AS4/3501-6 Unidirectional 
Continuous Fiber Reinforced Composite ( [v f =0.6) 


Elastic 

Constants 

(Gpa) 

Present 

results 

Analytical 
results (Ref. 
[10]) 

Analytical 
results (Ref. 
[11]) 

Experimental 
results (Ref. 
[12]) 

£, 

142.60 

142.9 

142.9 

142 

e 2 

9.51 

9.20 

9.79 

10.30 

G[i 

5.92 

5.50 

6.53 

7.60 

G 2 3 

3.05 

— 

3.01 

3.8 

^12 

0.25 

0.26 

0.26 

— 

^23 

0.35 

— 

0.42 

— 



Figure 2: Variation of Fiber Volume Fraction (Vy ) across the 

Thickness of Composite Famina for Different Values of 
Power Faw Exponent (V a f ) 
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Figure 4: Computed Effective Elastic Coefficients for Different Fiber 
Matrix Packs of RVE and the Corresponding Fitted Curves 
for Different Values of Sine Law Exponent 

(V a f ) (N = 3, vj = 0.6) 



Figure 5: Computed Effective Elastic Coefficients for Different 
Fiber Matrix Packs of RVE and the Corresponding Fitted 
Curves for Different Values of Sine Law Exponent 

(V a f ) (N = 3,v™ = 0.6) 


Table 3: Magnitudes of the Fitted Curve Parameters using Sine Law Variation 



a» 

al 

a2 

a3 

a4 

a5 

a6 

V f . 

0.15 

0.1 

0.15 

0.10 

0.15 

0.10 

0.15 

0.10 

0.15 

0.10 

0.15 

0.10 

0.15 

0.1 

Cll 

7.3000E+10 

7.0782E+10 

0 

0 

0 

-6.488 lE+11 

0 

0 

2.1000E+12 

2.1625E+12 

0.001 

0.003 

0 

-2.3497E+12 

C22 

1.2286E+10 

1.1976E+10 

0 

0 

0 

-5.8121E+10 

0 

0 

2.1058E+1 1 

2.0302E+11 

0 

0 

0 

-2.5166E+11 

C33 

1.2386E+10 

1.1966E+10 

0 

0 

0 

-5.7486E+10 

0 

0 

1.9477E+11 

1.9329E+11 

0 

0 

0 

-2.1587E+11 


2.8823E+09 

2.7630E+09 

0 

0 

0 

-1.3563E+10 

0 

0 

5.3374E+10 

4.9636E+10 

0 

0 

0 

-6.8721E+10 

C 55 

3.4244E+09 

3.3063E+09 

0 

0 

0 

-1.9055E+10 

0 

0 

6.3963E+10 

6.8381E+10 

0 

0 

0 

-9.0584E+10 

Cgg 

3.3928E+09 

3.3000E+09 

0 

0 

0 

-1.9000E+10 

0 

0 

6.2153E+10 

6.8000E+10 

0 

0 

0 

-9.0000E+10 

C 12 

5.0451E+09 

4.9705E+09 

0 

0 

0 

-1.6287E+10 

0 

0 

5.7801E+10 

5.7152E+10 

0 

0 

0 

-7.1683E+10 

Cl3 

5.0797E+09 

4.9831E+09 

0 

0 

0 

-1.6670E+10 

0 

0 

5.5097E+10 

6.2161E+10 

0 

0 

0 

-8.9593E+10 

C23 

5.3191E+09 

5.2632E+09 

0 

0 

0 

-1.8994E+10 

0 

0 

7.0034E+10 

6.2733E+10 

0 

0 

0 

-6.6224E+10 
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RESULTS AND DISCUSSIONS 


The effective elastic constants of the graded lamina are evaluated using the micro-mechanics finite element model 
of the corresponding RV. The present micro-mechanics finite element formulation and boundary conditions (Table 1) are 
verified considering uniform FVF (v f =0.6) in all fiber-matrix packs within the RV. The computed effective elastic 

constants of this continuous fiber-reinforced composite of uniform FVF are compared with those of an identical composite 
lamina analyzed by earlier researchers []. 

For two different values of the power-law exponent ( n r ), the magnitudes of a coefficient ( c if ) for different 

homogenized sub-volumes of RV are illustrated in Figure 2 and Figure 3. The variation of a coefficient along the thickness 
(Z) direction is mathematically modeled by a sine-function for the continuous variation of the fiber volume fraction as 
follows. 


i- \ m 

Vf(z) = Vf 


1 - 2 * 

h 


+ v°f sin(0)sin 


2 nz 


L sin(0) 


(16) 


6 = tan - :1 (h / (2vf )), / m = J(v' f n ) 2 + (h 2 /4) 

The Figure 1 shows the plot of different amplitudes versus fiber volume fraction, that with the increase in the 
amplitude of the sine function, FVF of the sub-volume decreases. Fig.2 and Fig.3 shows the computed magnitude of a 
coefficient for a homogenized sub-volume is plotted corresponding to the Z -coordinate of its (sub-volume) centre. For any 
of the two cases (n = 1 or 2), it may be observed from these figures that the nature of variation of a coefficient along the 

thickness (Z) direction is similar to that of the FVF (Fig. 1). The values of the fitted curves for different amplitudes of the 
sine function are given in Table 1. A lower value of FVF at the top and bottom layer is recommended for lower differences 
in the elastic properties at the top and bottom surfaces of the lamina. The value of amplitude and the number of sub- 
volumes considered are important in view of practical design of composites within a specified thickness. It can be observed 
that, increase in amplitude of the sine function yields a lower FVF at the top and bottom sub-volumes. Even though a 
continuous variation of material properties can be achieved in a composite laminate which is desirable in reducing inter 
laminar discontinuity, it may cause difficulties in the fabrication of the composite. 

CONCLUSIONS 

In the present work, a graded fiber-reinforced composite lamina is designed with an objective of reduced stress- 
discontinuity at the inter-laminar surfaces of traditional composite laminates. The representative volume (RV) of the 
graded lamina is composed of several micro-volumes stacked along the thickness direction. Every micro-volume has a 
continuous fiber and certain FVF. The FVF discontinuously varies among the different micro-volumes according to a 
simple power law. A continuum micro-mechanics finite element model of RV is derived and its effective graded elastic 
properties are predicted according to its (RV) actual constructional architecture. In this prediction, the classical 
homogenization treatment is implemented on every micro-volume without decoupling it from rest of the micro-volumes 
within the RV. The numerical results illustrate micro-volume-wise variation in the magnitude of every elastic coefficient. 
This variation for every elastic coefficient is also assumed as continuous variation according to a sine function. It is 
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observed that, a particular nature of variation of FVF yields the same nature of variation of every elastic coefficient across 
the thickness of RV. 
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